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We investigate the details of pattern formation and transitions between different modulated phases 
in ultra-thin Fe films on Cu(OOl). At high temperature, the transitions between the uniform satu- 
rated state, the bubble state and the striped state are completely reversible, while at low temperature 
the bubble phase is avoided. The observed non-equilibrium behavior can be qualitatively explained 
by considering the intrinsic energy barriers appearing in the system due to the competition between 
the short-ranged exchange and the long-ranged dipolar interactions. Our experiments suggest that 
the height of these energy barriers is related to the domain size and is therefore strongly temperature 
dependent. 
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I. INTRODUCTION 

Ultra-thin ferromagnetic iron films on the Cu(OOl)- 
surface are magnetized perpendicularly to the film 
plane 1 3 . Due to the competition between the short- 
ranged ferromagnetic exchange interaction and the long- 
ranged dipolar interaction, the magnetization breaks up 
into magnetic domains which are organized into more or 
less regular patterns^. The formation of patterns as 
a result from competing interactions on different length 
scales is observed in many physical and chemical systems 
ranging from type-I-superconductors in the mixed state 6 
or amphiphilic molecules at the air-water interface 7 to 
thermal convection 8 . Which pattern occurs depends on 
various parameters, in the case of magnetic films e. g. on 
the temperature T and the magnetic field H applied per- 
pendicularly to the film plane. In two dimensions, regular 
patterns include parallel stripes of alternating magneti- 
zation and circular bubble domains in a homogeneously 
magnetized background. Both types of patterns were ob- 
served in magnetic garnet films with a thickness of sev- 
eral micrometers^ - — . However, transitions between the 
stripe and bubble phase could only be induced if a small 
AC-field was applied to excite the system. In ultra-thin 
magnetic films metastable bubble patterns could be ob- 
served after applying a strong, almost in-plane, magnetic 
field pulse 14 . Also in the switching of magnetic multi- 
layer systems metastability and out-of-equilibrium states 
are the rule rather than the exception 15 . 

Recently, equilibrium transitions between stripe, bub- 
ble, and uniform patterns were observed at high enough 
temperature in ultra-thin Fe films on the Cu(OOl)- 
surface^ and the phase diagram in the T-i^-plane was 
experimentally determined. For convenience, the phase 
diagram is reproduced schematically in Fig. 1. In the 
present paper we investigate the details of the trans- 
formations between the different patterns by following 
paths of constant field or constant temperature in the 
phase diagram. While at high temperatures all pro- 
cesses are completely reversible, at low temperature we 
observe hysteretic behavior in both types of paths. More- 



over, we find that the bubble pattern is systematically 
avoided at low temperatures. These non-equilibrium as- 
pects result from energy barriers which intrinsically arise 
from the competition between the short- and long-ranged 
inter actions^ 7 -"—. They occur also under ideal conditions, 
i. e. in absence of pinning to structural defects. The en- 
ergy landscape of the system is explored in a ground- 
state continuum model, which qualitatively explains the 
observed behavior at all temperatures. 

The structure of this paper is as follows. In section II 
we discuss details of the sample preparation procedure, 
the experimental set up, and the general properties of 
our samples. In section III we introduce the important 
quantities that characterize a given sample. In section 
IV we image the transformations of the domain patterns 
along different paths in the T-i^-plane and discuss their 
(ir-) reversibility. In section V we present details of the 
transformations in the vicinity of the phase boundaries 
between different patterns. Finally, in section VI, we 
summarize the essential conclusions drawn from the ex- 
perimental observations. In appendices A and B, we 
present an extended report on the ground-state energy 
computations, some results of which are used to inter- 
pret experimental observations. 



II. SAMPLE PREPARATION AND 
EXPERIMENTAL SETUP 

Our samples are grown by molecular beam epitaxy 
at room temperature in ultra-high vacuum (UHV, base 
pressure below 1.5 • 10 -10 mbar) with a typical rate of 
0.1 to 0.2 atomic layers per minute. The substrate is 
a copper single crystal with the (OOl)-surface oriented 
to better than 0.1° and polished to < 0.03 fim rough- 
ness depth^l. The substrate is cleaned by several cycles 
of Ar-ion bombardment (argon pressure 2.5 • 10 -5 mbar, 
ion energy 1 keV) and annealing at 700 K for 50 min- 
utes. The chemical cleanliness and film thickness are 
determined using Auger electron spectroscopy, the crys- 
tallinity of the sample is verified by low-energy electron 



2 




M S (T) 




Lo(T) 






\ \ 
^-1.7- 







Fig. 4(a)-(e) 
— Fig. 4(f)-(k) 

uniform 
paramagnetic 

~-\g. 4(l)-(u) 



Nuc 



7~ c temperature 



FIG. 1. (color online) Schematic pattern phase diagram for 
ultra-thin Fe films with perpendicular magnetization accord- 
ing to Ref. 16. The gray double arrows illustrate the paths 
followed for the measurement of figures 2 to 4 as indicated. 
The vertical dashed lines mark characteristic temperatures 
Tnuc and Tc as discussed later in the text. The inset shows 
schematically the temperature dependencies of the important 
parameters Lo (blue-continuous curve) and Ms (red-dashed 
curve), see the text. 



diffraction. The film thickness d is measured with an ac- 
curacy of ±0.02 nominal atomic monolayers (ML). After 
preparation, the sample is transferred to the measure- 
ment chamber (base pressure below 3 • 10 -11 mbar) in 
the same UHV system. We measure the local magnetiza- 
tion of the sample with a lateral resolution of 50 nm using 
Scanning Electron Microscopy with Polarization Analy- 
sis of the secondary electrons (SEMPA)2i~— . By scan- 
ning the beam of a scanning electron microscope over 
the sample and analyzing the spin polarization of the 
secondary electrons generated at the electron beam spot, 
we obtain a spatial map of the local magnetization vec- 
tor. A magnetic field of up to 0.4 mT can be applied to 
the sample while measuring and the sample temperature 
can be varied from 20 to 400 K. The magnetic field can 
be controlled on the level of 1 fiT and the temperature 
is stabilized to less than 0.1 K. For electron-optical rea- 
sons, in our experiment the field is applied at an angle 
of 45° with respect to the film normal. However, within 
our experimental resolution we do not observe any in- 
plane component of the magnetization nor do we see an 
influence of the in-plane component of the field on the 
domain pattern when we rotate the sample about its ver- 
tical axis. We conclude therefore that for the very weak 
fields applied in the present experiments, only the out-of- 
plane component of the magnetic field affects the sample. 
Because the Fe films are very susceptible to weak mag- 
netic fields along their normal, it is important to elimi- 
nate the perpendicular component of any unwanted field 
(e. g. stray fields from the magnetic lenses of the electron 
microscope or the earth magnetic field) . This is achieved 
by tuning the applied field appropriately and zero-field is 
defined as the situation with no imbalance between up- 
and down magnetized domains. The magnetic field val- 



ues indicated throughout this paper refer to the effective 
perpendicular component of the field at the sample. In 
the images, the out-of-plane component of the magneti- 
zation is indicated by a gray scale with black and white 
corresponding to opposite signs of the magnetization. 



III. CHARACTERISTIC QUANTITIES 

In zero field, we observe a characteristic temperature 
Tc at which the contrast resulting from the magnetic do- 
mains disappears in the SEMPA images 2 ^. This sample- 
specific Tc sets the temperature scale for each individ- 
ual film. By introducing a reduced temperature T/Tc, 
the temperature dependence of the domain patterns can 
be compared between different samples with slightly dif- 
ferent film thickness (ranging from 1.9 to 2.3 ML) and 
different Tc (ranging from 320 to above 360 K). As a 
function of the film thickness, Tc reaches a maximum for 
2.15 ML, decreasing rapidly for thinner films and more 
slowly for thicker films 25 . In zero field, the domain pat- 
tern consists of parallel stripes of width Lo with alternat- 
ing magnetization ±Ms as expected 1 . Throughout the 
temperature range considered in this work, we observe a 
strong decrease of L$ with increasing temperature 2 ^ 2 ^, 
see the inset in Fig. 1. As discussed in appendix A, this 
stripe width in zero field, L , determines the domain size 
at all values of the applied field and is therefore an im- 
portant quantity characterizing the sample at any given 
temperature. If a magnetic field is applied, the domains 
with the appropriate sign of the magnetization are ener- 
getically favored and their area grows at the expense of 
the domains carrying the opposite sign. The pattern ac- 
quires an asymmetry A= (f t - f±)/ (/ t + where 
is the area occupied by up (down) magnetized domains. 
The zero-field pattern with alternating stripes of equal 
width has A = and the up or down magnetized satu- 
rated states in the limit of high fields have A = ±1. For 
intermediate values of A at intermediate magnetic fields 
one expects that the minority, down-magnetized domains 
form a hexagonal lattice of (almost) circular (bubble) 
domains in a homogeneous up-magnetized background 
rather than a striped pattern 2 ^. Because the weak mag- 
netic fields considered in this study only change the do- 
main configuration but not the local magnetization in- 
side the domains, the latter is equivalent to the satu- 
ration magnetization, Ms- The total magnetization of 
the sample for any value of A is then given by the prod- 
uct A- Ms- Over most of the temperature range Ms is 
weakly temperature dependent, except close to Tc where 
it drops sharply, see the inset in Fig. 1. The magnetic 
field scale for a given sample at a given temperature is 
determined by its critical magnetic field He, defined as 
the field above which the uniform state \A\ = 1 has the 
lowest energy. It turns out that the asymmetry A is pro- 
portional to the reduced field H/Hc over most of the 
magnetic field ranged. Moreover, He can be expressed 
through the saturation magnetization Ms and the stripe 
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width in zero field L , He cx M^/Lq, as discussed in 
appendix A, Eq. (|A7|) . Together with the experimental 
curves for M S (T) and L (T)^^^, this relation allows 
us to trace the phase boundary separating the uniform 
and the bubble states in Fig. 1, 

Inserting typical values for d = 2.0 ML = 0.36 nm, 
HqMs = 2 T and L$ = 1 jam corresponding to our ex- 
perimental system a few K below Tc, we obtain ji^Hc = 
390 //T, in good agreement with the typical saturation 
fields observed in the experiments, see the next sections 
and Ref. [l6|. 



IV. DOMAIN PATTERNS VS. FIELD AND 
TEMPERATURE 

A. Constant temperature 

The most obvious way to observe a transition from 
stripe to bubble domains is to apply a magnetic field at 
constant temperature. For T close to Tc, this process is 
shown in Fig. 2. We start from a striped pattern at al- 
most zero field, image 2(a), with white and black regions 
almost balanced. In image 2(b), the external field results 
mainly in a compression of the (white) minority stripes. 
In this step splitting of selected domains has taken place 
and a few domains have acquired a bubble-like shape. 
This trend is continued in image 2(c), where the numbers 
of stripe-like and bubble-like domains are roughly equal. 
In image 2(d) virtually all stripe domains have broken up 
and the domain pattern consists of a disordered array of 
white bubbles immersed in a black background. Upon in- 
creasing the magnetic field further, the bubble diameter 
decreases but the bubble density remains approximately 
constant. The transition to saturation is not observed in 
this case because imaging is not feasible with our experi- 
mental setup in the fields required to saturate the sample 
at this temperature. Images 2(f) to 2 (j) show the domain 
transformation in decreasing magnetic field, still at high 
temperature. The barely visible bubble domains of im- 
age 2(f) partially merge, leading to randomly distributed, 
slightly elongated domains in image 2(g). These merge 
further and in image 2(h) only a few bubble domains are 
left and the pattern clearly has acquired a stripe charac- 
ter. Finally, in images 2(i) and 2(j) the stripe pattern is 
restored completely. 

Note that the periodicity of the domain pattern, i. e. 
the typical distance between neighboring bubbles or 
stripes with the same sign of the magnetization, is es- 
sentially unaffected by the magnetic field. This obser- 
vation is in agreement with ground state energy calcula- 
tions, see Fig. 8 and the appendix A. Since the period- 
icity of the pattern does not change, the transition from 
black-and-white stripes to white bubbles on a homoge- 
neous black background implies that the total number 




FIG. 2. Bubble formation upon changing the magnetic field 
at constant temperature. The absolute value of the magnetic 
field for each image is indicated. Images (a)— (e): increas- 
ing magnetic field at T/T c =0.994 (T=350 K, d=2.00 ML), 
(f)-(j): decreasing magnetic field at T/T c =0.996 (T=339 K, 
<i=1.93 ML). All image sizes are 45 /mi by 45 /xm, the inset 
illustrates the path followed in the phase diagram. 

of white domains increases and splitting of the white do- 
mains must occur. Alternatively, we could also state that 
the black stripes have to merge to form the homogeneous 
background. Merging and splitting are therefore com- 
pletely equivalent - which process occurs depends only 
on whether one describes the black or the white domains. 
For the sake of simplicity from here on we concentrate on 
describing the white domains which constitute the minor- 
ity domains for the direction of the magnetic field applied 
here. As discussed in more detail in the appendix B, split- 
ting of domains implies the crossing of an energy barrier 
because the energy cost associated with temporarily in- 
creasing the length of the domain walls is not completely 
compensated by the gain in dipolar energy. The exis- 
tence of this energy barrier is an intrinsic property of the 
pattern forming system and has nothing to do with the 
extrinsic energy barriers responsible for the pinning of the 
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FIG. 3. Domain patterns upon changing the magnetic field at 
constant T=317 K (T/T c =0.970). The absolute value of the 
magnetic field is indicated for each image. Images (a)— (e): 
increasing field, (f)— (j): decreasing field. Note that images 
(f)-(j) have been recorded before images (a)-(e) with (j) and 
(a) being the same image. All image sizes are 45 /mi by 45 /mi, 
the inset illustrates the path followed in the phase diagram. 



domain walls at structural defects such as atomic steps 
on the substrate or the sample border. The fact that the 
sequences 2(a)-(e) and 2(f)-(j) are reversible indicates 
that for comparatively high temperatures merging and 
splitting readily occur, and the system reaches its equi- 
librium stripe, bubble or uniform state in agreement with 
the phase diagram of Ref. [l6l (Fig. 1). Moreover, pinning 
at minor structural defects seems to play no role in this 
case. 

At lower temperatures (Fig. 3) domain splitting is more 
difficult because the system has less thermal energy to 
overcome the energy barrier and additionally, the energy 
barrier itself increases as discussed in the appendix B. 
There we suggest that the height of the energy barrier 
which impedes domain splitting or merging is propor- 
tional to the characteristic domain size Lo- This implies 
that the energy barrier also has the same temperature 



dependence as Lq(T), meaning that the barrier strongly 
increases for decreasing T. Temperature-dependent en- 
ergy barriers have been suggested in a frustration-based 
approach to supercooled liquids in the context of glass- 
forming systems^. The fact that splitting is forbidden 
implies that the total number of domains in the sample 
can only increase by nucleation of new domains. There- 
fore the response of the system to the applied field is 
fundamentally different from the equilibrium situation of 
Fig. 2. Applying a weak field to the stripe pattern of 
image 3(a) results in a lateral compression of the stripe 
domains, see image 3(b), in agreement with the expecta- 
tions from the ground state calculations and the equilib- 
rium situation. For larger fields, the system can not reach 
the equilibrium bubble state because this would require 
domain splitting and in consequence the system is locked 
in a metastable stripe pattern, image 3(c). Nevertheless, 
the stripes successively disappear, image 3(d), finally giv- 
ing way to the uniform saturated state, image 3(e), which 
is unavoidable at large enough fields. Because a large en- 
ergy barrier prohibits the lateral collapse of the stripes, 
see the appendix B, a stripe can only disappear by con- 
tracting along its length. This is not possible as long as 
the ends of the stripe are pinned at some defect or the 
sample border. At a high enough field, however, we ex- 
pect that one end of the stripe will unpin, followed by 
a rapid contraction of the stripe until it is reduced to a 
single circular bubble domain. Such residual bubble do- 
mains can be observed in images 3(b), 3(c), and 3(d). As 
it is pointed out in the appendix, the energy barrier that 
stabilizes the isolated bubble domain disappears above a 
critical field i^coiiapse = 1.35ifc? where He is the equi- 
librium saturation field introduced in section III. There- 
fore, for H > i^coiiapse the bubble inevitably collapses, 
without crossing an intrinsic energy barrier. In conse- 
quence, this process is allowed also at low temperature, 
at least as long as domain wall motion in general is pos- 
sible. According to this discussion, the bubble phase is 
systematically avoided at low temperature, even in the 
ideal system, in absence of pinning at structural defects. 
The fact that the images in Fig. 3 have been recorded 
only 9 K below Xc, illustrates the difficulty of observing 
the equilibrium bubble pattern at constant temperature 
in this system. 

When reducing the magnetic field from the saturated 
state at low temperature, image 3(f), only a small num- 
ber of domains nucleate, image 3(g), probably at some 
structural defects that are not resolved in the image 28 . 
Nucleation of a domain - in contrast to its collapse - is 
always hindered by an energy barrier, see the appendix 
B, although a local modification of the film properties 
at a structural defect may lower this barrier and pro- 
mote domain nucleationi 2 -. Note that image 3(g) has 
been recorded at the same magnetic field as image 3(b). 
The initial seed domains in image 3(g) expand along one 
direction, image 3(h), and upon lowering the field further, 
more stripe domains enter the field of view, image 3(i), 
until the weight of up and down magnetized domains is 
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essentially equilibrated in image 3(j). 

In summary, while at high temperature, in close prox- 
imity of Tc, the observed domain pattern closely follows 
the equilibrium phase diagram of Ref. 16, already a few 
Kelvins below Tc the domain pattern is dominated by 
metastability. The transition from the striped to the 
saturated state in increasing field does not proceed via 
the bubble phase but by direct collapse of the stripe 
domains 15 . 



B. Constant applied field 

As is evident from the phase diagram of Fig. 1, tran- 
sitions between the different phases can also be induced 
by varying the temperature while keeping the applied 
field constant. In Fig. 4 we present the domain pat- 
terns observed by following such paths of constant field, 
as indicated in the inset. The images of Fig. 4 are ob- 
tained by scanning the electron beam over the sample 
while the temperature is varied at a fixed rate. With 
a cooling/heating rate of ^0.5 K/min and acquisition 
of each image taking 400 seconds, each one of the im- 
ages 4(a) to 4(p) comprises a temperature variation of 
3.2 K. A different position in the image is therefore equiv- 
alent to a different sample temperature. The temper- 
ature scale is indicated in terms of the reduced tem- 
perature by the vertical lines in Fig. 4 and the images 
are placed accordingly. Note that images 4(q) to 4(u) 
have been acquired at a slightly lower heating rate of 
0.4 K/min and have been scaled by a factor of 0.8 along 
the horizontal axis in order to fit on the same tempera- 
ture scale. 

For high magnetic fields [241 /xT, images 4(a)-4(e)], 
no stripes are observed, in agreement with the phase di- 
agram of Fig. 1. When cooling from above Tc, Fig. 4(a), 
bubble domains appear in image 4(b) whose density de- 
creases monotonously with decreasing temperature un- 
til a homogeneous, saturated state is reached, Fig. 4(c). 
Since for this field the transitions occur at high tempera- 
ture, the process is completely reversible, as can be seen 
in images 4(d) and 4(e) recorded upon heating. We can 
therefore be sure that the system is at equilibrium, as in 
Fig. 2. 

For intermediate values of the magnetic field, im- 
ages 4(f)-4(h), the sequence of phases upon cooling is 
uniform paramagnetic — » stripes — > bubbles — » uniform 
saturated. The transition from stripes to bubbles occurs 
at a high temperature of about 0.99 Tc, see images 4(f) 
and 4(g), at which domain splitting occurs and the pro- 
cess is equivalent to the observations of Fig. 2. Upon 
decreasing the temperature further, the density of the 
domains decreases by successive collapse of the bubbles, 
until the saturated state is reached in image 4(h), at a 
clearly lower temperature than for the high field in im- 
age 4(c). Remember that collapse of the bubble domains 
is allowed also at low temperatures. Upon heating, a 
different behavior is observed. The sample remains satu- 



rated until a sample-specific temperature Tn uc is reached 
and suddenly, from one scan-line to the next one, stripe 
domains appear on the sample, image 4(i). As pointed 
out before, nucleation always implies the crossing of an 
energy barrier, for which a certain thermal activation is 
required. In the appendix B we compute the height of 
this nucleation barrier Eb for a single bubble domain 
in the ground state. Analogously to the discussion of 
Hc(T), see Eq. (pQ) in section III, we may translate the 
ground-state result of Eq. (|B7j) to finite temperatures 
by using the experimental results for Ms(T) and Lq(T). 
For H = He corresponding to the equilibrium transition 
from the uniform to the bubble state we obtain 



E B (T) = 0.1904/i Ms(T) 2 d 2 To(T) . (2) 



Eq. (j2j) suggests that the barrier height Eb is propor- 
tional to To. Using the same values as before, d — 
0.36 nm, n^Ms = 2T and To = 1 /im, in (|2]), we obtain 
Eb/^b = 5680 K, where ks is the Boltzmann constant. 
At an absolute temperature of 330 K, we obtain then a 
probability exp (— Eb/^bT) = 3.3 • 10 -8 , which gives a 
reasonable nucleation rate if one assumes an attempt fre- 
quency of O (10 9 - 10 12 ) Hz 2 ^. Note that according to 
([2]), the main temperature dependence of the nucleation 
probability is not due to the varying absolute tempera- 
ture T, but rather due to the strong temperature depen- 
dence of Eb(T) via To(T). The nucleation temperature 
Tnuc is therefore determined by the reduced temperature 
T /Tc at which To(T) (and with it Eb) drops below some 
threshold value. The absolute temperature plays only a 
minor role. 

The stripes in image 4(i) decay to the equilibrium bub- 
ble state, Fig. 4(j), as soon as the temperature is high 
enough to allow for the splitting of domains. At higher 
temperature the equilibrium transition from bubbles to 
stripes occurs, and the stripes in turn disappear as Tc 
is reached, Fig. 4(k). The transition from saturation to 
bubbles via transient stripe-like domains is discussed in 
more detail in the next section in the context of figure 5. 

At lower fields, the equilibrium phase sequence is again 
observed upon cooling. Because the temperature inter- 
val for each phase is wider and the domain size at lower 
temperature is larger, the phases can be identified more 
clearly in the images 4(l)-4(p). Upon heating, the cross- 
ing of the nucleation temperature Tn uc can be identi- 
fied clearly, although a few domains are visible already 
at lower temperature in images 4(q) and 4(r). In this 
case Tnuc lies above the bubble-stripe transition tem- 
perature, compare images 4(n) and 4(r). Therefore no 
bubble phase is observed upon heating in low fields, the 
nucleated stripe pattern is stable and persists upon fur- 
ther heating up to Tc- 
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FIG. 4. (color online). Bubble formation upon changing the temperature in constant field. The images (a)-(p) were acquired 
at a constant cooling/heating rate of 0.5 K/minute, leading to a continuous temperature variation of 3.2 K within each image. 
Image (a)-(c): cooling in 241 fiT, (d), (e): heating in 241 /xT, (f)-(h): cooling in 195 fiT, (i)-(k): heating in 195 /iT, 
(1)— (p): cooling in 146 /iT, (q)— (u): heating in 146 /iT at a rate of 0.4 K/min. The images (q)-(u) have been scaled by a 
factor of 0.8 along the horizontal direction to fit on the same temperature scale. All image sizes are 45 /im by 45 /im. The film 
thickness is <i=2.22 ML, all images were measured on the same sample. The vertical lines indicate the reduced temperature 
scale and the inset illustrates the paths followed in the phase diagram. 



V. DETAILS OF THE PATTERN 
TRANSFORMATIONS 

At high applied fields and rising temperature the tran- 
sition from saturation to the bubble pattern proceeds by 



direct nucleation of circular domains. At low fields the 
bubble phase is avoided completely and nucleation di- 
rectly leads to a striped pattern. The transition from 
the saturated state to the bubble state via nucleation of 
a transient stripe pattern as observed in Fig. 4(i)-4(j) 
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FIG. 5. (color online). Transition from the saturated state to 
the bubble state via a transient, string- like domain pattern in 
a constant magnetic field of 156 /jlT. During the acquisition 
of image (a) the sample is heated from 332.5 K to 335.5 K. 
In image (b) the temperature settles to 336.5 K and is then 
held constant for the remaining images. All image sizes are 
45 fim by 45 /mi, the film thickness is <i=2.00 ML. The slow 
scan direction is from top to bottom and the acquisition time 
is 3.5 minutes for each image. 



is presented in more detail in Fig. 5. An initially sat- 
urated sample is heated in a constant magnetic field of 
156 jaT with a constant rate of +1 K per minute. When 
the nucleation temperature marked by the red dotted 
line in image 5(a) is crossed, suddenly reversed white 
domains appear on the sample. This process can be un- 
derstood as follows. The uniformly magnetized state be- 
comes metastable when the equilibrium transition line 
from bubbles to saturation is crossed upon heating. How- 
ever, nucleation is prohibited until the associated energy 
barrier has sufficiently decreased, as discussed in the pre- 
vious section and the appendix B, and at some temper- 
ature Tnuo isolated bubble domains nucleate. Since the 
expansion of an existing isolated domain is not hindered 
by an energy barrier, these domains grow rapidly (instan- 
taneously, on the time scale of our measurement) until 
the equilibrium asymmetry A corresponding to the ap- 
plied field and temperature is reached. This process is in 
agreement with the observations by Cape and Lehman^ 
in thick garnet films. At this point the temperature ramp 
is halted and T is kept constant during the remaining im- 
ages. As can be seen from the contiguous domain marked 
yellow in image 5(b), the initial domains extend over 
large areas. However, in contrast to thick garnet films, in 
our case thermal fluctuations are strong enough to break 
up the string-like domains and increase the total num- 
ber of domains. This process can be observed already 
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FIG. 6. (color online). The sequence of SEMPA images shows 
the domain pattern transformation on heating a bubble do- 
main pattern from 300 K up to 350 K. The magnetic field is 
22 fiT. Image size is 182 fim by 182 fim for all images and the 
film thickness is <i=2.16 ML. 



in image 5(c) where the many colored domain segments 
originate from the single domain marked in image 5(b) 
and is more pronounced in images 5(d), 5(e), and 5(f). 
Finally, in images 5(g) and 5(h) only bubble domains of 
approximately circular shape are left. 

The transformation from bubbles to stripes in a con- 
stant field of 22 /iT is investigated in Fig. 6. To help 
following the pattern transformation, one representative 
domain is highlighted in image 6(a) and the same do- 
main is marked yellow in all images. Images 6(a) and 
6(b) show that at low temperature the domain pattern 
is mostly unaffected by the temperature change of 30 K. 
Both images show a pattern of essentially randomly ar- 
ranged round domains. In image 6(b) the size distri- 
bution is narrower and the domains have a more circu- 
lar shape, indicating that the pattern is closer to equi- 
librium. From image 6(b) to 6(c) the bubble domains 
elongate and no preferential direction can be observed 
for this elongation. In image 6(d) the aspect ratio of 
the domains increases further and in image 6(e) the do- 
mains start to preferentially align along the same direc- 
tion. Note that in spite of the high absolute tempera- 
ture, no domain merging or splitting occurs due to the 
low reduced temperature (Tc > 360 K) and associated 
large domain size. Therefore the total number of domains 
remains constant, the transformation from the bubbles 
to the stripes is achieved only by the growth in length 
and shrinking in width of the domains. This holds true 
also for image 6(f), in which the order of the stripe pat- 
tern increases further. The transformation from bubbles 
to stripes in increasing temperature at constant field is 
driven by two effects. Because the equilibrium periodic- 
ity of the pattern Lo decreases, the width of the white 
domains tends to decrease. However, because Hc(T) in- 
creases and A oc H/Hc , the equilibrium asymmetry A 
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FIG. 7. (color online). The sequence of SEMPA images shows 
the domain pattern transformation on cooling from above the 
Curie temperature Tc to saturation. The magnetic field is 
95 /iT and all image sizes are 45 fim. by 45 /mi. The film 
thickness is d=1.91 ML. 



decreases and therefore the area occupied by the white 
domains tends to increase. The decrease in width and 
simultaneous increase in area favors an elongated shape 
of the domains and thus a transformation from bubbles 
to stripes with increasing temperature, without the need 
for domains to merge. 

A detailed measurement of the domain pattern evolu- 
tion upon cooling from above Tc in a constant applied 
field is shown in Fig. 7. The sample is cooled in discrete 
steps of -1 K and after the temperature is allowed to set- 
tle during 3 minutes, the acquisition of the next image is 
started. All images show the same field of view. The im- 
age 7(a) taken at 330 K w Tc is essentially contrastless 
and shows only small regions of stripe domains marked in 
yellow and blue. In image 7(b) a mostly striped pattern 
can be observed, intermixed with some bubble domains. 
The next image, 7(c), shows a mixture of elongated do- 
mains and approximately circular bubbles with a larger 
average domain size than in image7(b), in agreement with 
the larger Lq at lower T. At a temperature of 327 K, in 
image 7(d), except for a few all domains have a circular 
shape and in image 7(e) only circular domains are left. As 
it is obvious from the image, the bubbles are not arranged 
on a hexagonal lattice as one might expect but they are 
rather disordered as was observed already in Fig. 5(d) 
and 5(e). Since the fast scanning direction of the micro- 
scope is from left to right, the fact that domains are mov- 



ing during scanning manifests as horizontal streaks in the 
image. Two domains showing such streaks are marked 
yellow in image 7(e). This mobility is also observed in 
image 7(f), for example in the yellow colored domains. 
From images 7(e) to 7(f) and 7(g) the bubble density de- 
creases with decreasing temperature. In image 7(h) only 
two reversed domains are left and finally, in image 7(i), 
the entire field of view is homogeneously magnetized. As 
long as domain merger, splitting, and nucleation are not 
possible, the total number of domains in a sample can 
only decrease with time, since only domain deformation 
or collapse are allowed. For decreasing temperature, the 
increase in the reduced field H/Hc(T) tends to decrease 
the area occupied by the white domains. This decrease 
of the area proceeds by contraction of the stripes along 
their length. Simultaneously, the increase in Lq(T) favors 
an increase of the domain width and domain separation. 
Therefore, the temperature dependencies of both, Lq(T) 
and H/Hc(T) favor a transition from elongated to cir- 
cular domains with decreasing T, analogously to the case 
of increasing temperature discussed in Fig. 6. 



VI. CONCLUSIONS 

Our experiments show that in two-dimensional, ultra- 
thin ferromagnetic films with perpendicular anisotropy 
thermal fluctuations alone are sufficient to promote the 
transitions between stripe, bubble, and saturated phases 
in a narrow temperature interval below Tc. However, 
in general the actual pattern realized at a given mag- 
netic field and temperature depends not only on the equi- 
librium pattern as deduced from the phase diagram in 
Fig. 1, but also on the history of the sample and the spe- 
cific path followed in the T-i^-plane. At constant tem- 
perature, the characteristic size of the domain pattern 
is constant and a transition between the stripe and the 
bubble pattern necessarily involves merger or splitting of 
the domains. These processes, as well as nucleation, are 
inhibited by energy barriers that are inherent to the ideal 
system and are a direct result of the competition between 
the short- and long-ranged interactions. Translating the 
results of our ground-state computations to finite T sug- 
gests that the heights of these energy barriers are pro- 
portional to the domain size Lo and do therefore strongly 
depend on temperature, in the same way as does Lq(T). 
We conjecture that for decreasing temperatures domain 
merger, splitting and nucleation are suppressed not only 
because there is less thermal energy available, but mostly 
because the associated energy barriers increase strongly. 
In contrast, on a path of constant applied field H and, 
for instance, decreasing temperature, both the increase of 
Lq(T) and the increase of H/Hc(T) promote a transition 
from a stripe to a bubble pattern. Moreover, in this case 
the transformation is achieved by changing the shape of 
the existing domains only, no domain splitting is required 
and no energy barrier needs to be crossed. As a summa- 
rizing rule we may state that below a certain threshold 
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temperature the number of domains in the sample can 
only stay constant or decrease with time. While a pure 
rearrangement of the domains may accommodate a tran- 
sition from stripes to bubbles and vice versa for changing 
temperature, this rule forbids the transition from stripes 
to bubbles at constant low temperature. The pinning to 
structural defect in the sample and on the substrate plays 
no important role in the samples and in the temperature 
range considered in this work. It is therefore tempting to 
speculate that metastability and the associated disorder 
in the domain pattern are indeed self-generated by the 
long-range interaction^"— and entropy. 
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Appendix A: Ground-state phase diagram 

We compute the total energy of the domain patterns 
in different situations by minimizing the total energy of 
a slab with thickness d along the ^-direction and lateral 
extension A ^> d in the x-?/-plane. In a continuum ap- 
proach, we may represent the configuration of magnetic 
domains by a two- valued scalar field m(p) = ±1 if we 
assume that the magnetization points only along ±z and 
is homogeneous along the small thickness d. With this 
we have assumed that the domain walls are sharp and 
the local magnetization at a position p in the a>?/-plane 
is given by Msm(p)e z . The Hamiltonian for this model 
reads 

2Jdl wall + \JJ m(p)V d (\p-ff\)m((f)d 2 pd 2 p' 

-hd J m(p)d 2 p. (Al) 

The first term, corresponding to the exchange interaction 
with coupling strength J per atom, increases the energy 
stored in the domain walls by an amount proportional 
to the film thickness d and the total length of the do- 
main walls in the sample l wa ii- All lengths are measured 
in units of the lattice constant. The second term repre- 
sents the dipolar energy (coupling constant A = ^M|) 
as computed from the interaction between the magnetic 
surface charges at the top and bottom surfaces of the slab 
following Refs. I9I-H1I with 



V d {\P-f?\) 



1 



1 



\p-?\ y/\?-?\ 2 + # 



(A2) 



as given in more detail in Ref. [3l|. The last term in 
(jAl|) describes the interaction with the external field, 
ft = poHMs, and is equivalent to —dA 2 poHMsA(H), 



with A = -£2 f m(p)d 2 p being the asymmetry as in- 
troduced in section III. We consider the following do- 
main configurations: i) parallel stripes along the ?/-axis, 
with alternating magnetization ra = ±Ms and alternat- 
ing width L ± i, such that the periodicity of pattern 
is Ls = 2L and the width of the narrower stripes is 
w = L — 5. ii) a hexagonal array of reversed (ra = — Ms) 
bubbles of radius R and center spacing ona homoge- 
neous background (ra = +M5). hi) the saturated state 
ra = For a given set of the parameters J, <i, A, 

and ft > we minimize (|A1|) with respect to L and w or 
L b and i?, respectively, and compare the energies of the 
optimized stripe, bubble, and uniform patterns. 

For ultra-thin films we expect d <C (Ls,w,Lb,R)- 
In this limit, the energy per sample volume of the 
stripe configuration has been calculated by Kashuba and 
Pokrovsky 5 . Using our model we obtain an essentially 
equivalent analytical expression^, 



es 



2J 
L 



4Xd 



In ( — : sin C^^\ 



7rd 



\2L J 



t> (A3) 



except for the constant | coming from the slightly dif- 
ferent treatment of the dipolar energy and where we 
have subtracted the energy of the homogeneous, satu- 
rated state. For any ft, the w and L minimizing (jA3|) 
and the minimized energy can be found analytically. For 
ft = 0, w is equal to L and by minimizing (jA3|) w. r. t. L 
we obtain the stripe width in zero field, 



T ^ ( 2J 

Lo = d 2 e ^{4Xd 



(A4) 



For the bubbles, the energy per unit volume can be ex- 
pressed in a series expansion 31 



eB 



4ttR 
V3L 2 



2 J + Rh 



4X d\n(^=) 



\d^J 



00 / R \ 
+67rAd^S fc (— J 



A;=0 



3+2/c 



(A5) 



where the first terms correspond to the energy of a single 
bubble domain as given by Thiele 32 and the sum accounts 
for the interaction with the other bubbles in the sample. 
Note that the coefficients Sk are independent of both, R 
and Lb, and need to be computed only once. Because 
yt < 4, the series converges fast and for the precision 

B 

required here it is sufficient to consider only the first 10 
terms of the series. For high magnetic fields Lb ^> R 
and in this limit the series vanishes. By minimizing the 
remaining terms with respect to R and equating the min- 
imum energy to zero, we can find analytically the critical 
field above which the uniform state has the lowest energy, 



he = 32A exp 



2 J 3 
~4\d 2 



I6 71 " , d s 
-2- A — . A6 

e z L 



If we substitute the definitions of ft and A back into (|A6ft , 
we find a simple scaling expression relating the physical 
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FIG. 8. (color online), (a) Energies of the optimal stripe 
and bubble patterns (red-dashed and blue-continues curves, 
respectively) as a function of the applied field h, left scale. 
Black solid line: energy difference es — es, right scale. The 
vertical lines indicate the cross-over fields hs^B = 0.545/ic 
(short dashed) and hs^u = 0.924/ic (long dashed), (b) Size 
of the optimal domain patterns vs. applied field in units of 
Lq. Upper curves: pattern periods Ls and Lb (red-dashed 
and blue continuous, respectively). Lower curves: narrow 
stripe width w and bubble radius R (red-dashed and blue 
continuous, respectively). Vertical lines as in (a). 



critical field He to the characteristic quantities Ms and 
L . 

Hc = ^ r = - 2 M S ±. (AT) 
Ho Ms e 2 L 

In contrast to the stripes, for < h < he the energy of 
the bubbles (jA5|) can only be evaluated numerically. 

For the following analysis we fix the parameters J = 1 , 
d = 1, and A = 0.1, giving L = 141.4 and h c = 0.00481. 
Fig. 8(a) shows the energies of the stripe and bubble 
patterns (red dashed and blue continuous lines, respec- 
tively, left scale) obtained by minimizing the expressions 
(|A3|) and (|A5|) w. r. t. (L, w) and (L#, R), respec- 
tively. As it is clear from the graph, both energies in- 
crease monotonously for increasing fields and they differ 
at most by a few percent. For better clarity, the dif- 
ference between the two energies (e# — es, black solid 
line, right scale) is also plotted. The two energies be- 
come equal at a critical field hs^B — 0.545/ic (vertical 
short-dashed line), below which the stripes are favored 




Stripe width w/L Bubble radius R/L 



FIG. 9. (color online), (a) Energy per unit length of a single 
reversed stripe domain as a function of its width w for differ- 
ent values of the applied field /i, which is indicated for each 
curve in units of A. (b) Energy of a single reversed bubble 
domain as a function of its radius R for different values of h 
in units of A, as indicated for each curve. 



and above which the bubbles have a lower energy. We 
find numerically that the value of hs^B/hc = 0.545 
is independent of the specific choice of J, and A, 
as long as d/Lo <C 1. At a second characteristic field 

2 

hs^u = \hc = 0.924/ic (vertical long-dashed line), the 
energy of the stripes becomes equal to the energy of the 
saturated state, while the bubble pattern still provides 
an energetic advantage. The value of hs^u can be com- 
puted analytically by taking the limit L — >> oo in (|A3[) . 
see Eq. (jBl|) below. Hence, all characteristic fields are 
proportional to h c and therefore obey the same relation 
H ~ Msd/Lo, see Eq. ((All) and Ref - M 

In Fig. 8(b) we plot the characteristic lengths of the op- 
timized patterns. The two upper curves indicate the peri- 
odicities of the patterns Ls and L b (red dashed and blue 
continuous, respectively), while the two lower curves rep- 
resent the size of the minority domains, w and R. Note 
that the latter are almost equal throughout the entire 
field range and they remain finite at their respective crit- 
ical fields hs^u and he- The transitions to the saturated 
state -r 1 > and 7^ y do not occur through a van- 

Ls Lb 

ishing numerator, but by a diverging denominator. Nev- 
ertheless is the variation in the periodicity rather small 
throughout most of the field range, except close to the 
critical fields hs^u and he- 



Appendix B: Energy barriers and metastability 

Now we want to have a closer look at the mestastability 
of isolated domains around and above he as well as the 
energy barrier impeding domain merger and splitting in 
a domain pattern. 

Fig. 9(a) shows the energy (per unit length) of a single 
reversed stripe domain in an otherwise homogeneous in- 
finite film as a function of its width w, for different values 
of the applied field h. The value of h is indicated in units 
of A for each curve. The expression for the energy per 
unit length is easily obtained by multiplying Eq. (|A3|) by 
the periodicity Ls = 2L and the sample thickness d and 
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then taking the limit L oo, 



height E B 



stripe 



4 Jd - SXd 2 



In 



2hwd. (Bl) 



For h > the stripe energy has a shallow minimum at 
some value w* . For w —> Eq. (|B1[) contains an un- 
physical divergence and is no longer valid. Instead, the 
stripe reduces to two adjacent domain walls with a total 
energy corresponding to the associated cost in exchange 
energy of AJd per unit stripe length. This provides a 
large energy barrier that is independent of h and there- 
fore the isolated stripe domain remains metastable, even 
for values of h considerably exceeding the equilibrium 
saturation field he- Since the stripe domain can not col- 
lapse laterally, i. e. by letting w —> 0, eventually it will 
unpin at its ends or split at some point and then contract 
along its length. 

The expression for the energy ^bubble of an isolated 
bubble domain 32 with radius R and magnetization -\-Ms 
in a homogeneous background of magnetization — Ms is 
recovered from Eq. (|A5|) by multiplying the latter by the 
volume of the unit cell of the bubble lattice, d^-L 2 B , and 
then letting Lb — » oo which corresponds to neglecting 
the sum, 



^bubble = 2nRd 



2J-4Adln (J^Q +Rh 



(B2) 



This energy is plotted in Fig. 9(b) for d = 1, A = 0.1, J = 
1 and some values of the applied field h as indicated in 
units of A for each curve. For h small but > 0, the 
energy of the bubble domain has two local minima, one 
at R = and one at some R*. For h > he, the state 
with the bubble domain (R > 0) becomes metastable, 
but it is separated form the uniform state (R = 0) by an 
energy barrier. However, above a critical field /^collapse = 
| foes the local minimum at R* disappears and the bubble 
domain inevitably collapses to R = 0. This behavior 
is fundamentally different from the stripe domains and 
we expect that, as long as the domain walls are free to 
move in the sample, there will be no bubble domains 
for h > ^collapse? while isolated stripe domains whose 
ends are pinned at some defect or the sample border may 
persist up to considerably higher fields. 

The reverse process moves the system from the state 
R = to R = R* and corresponds to the nucleation 
of a single bubble domain. For any h this process is 
hindered by an intrinsic energy barrier as can be seen 
in Fig. 9(b). The local maximum of the energy defining 
the height of the nucleation barrier is found by setting 
dEbubbie/dR = 0. By using this condition and Eq. (jA6|) 
in (|B2|) we obtain the following expression for the barrier 



E B = 87r\d 2 L 



R-n 



Ln 



h A-k R mi 
h c e 2 L 



(B3) 



where i? max is the position of the local maximum of the 
energy, corresponding to the barrier height. The value of 
the ratio R max /Ln is obtained by using the definition of 
Lo, Eq. (|A4[) . in the condition <9L?bubbie/<9i? = 0, leading 
to 



ln ( 



47rR n 



h 2 47rR m 
h c e 2 L 



0. 



(B4) 



Note that the R max /Lo is only a function ofh/hc- Hence, 
the energy barrier impeding domain nucleation follows 
the scaling relation 



E B a\d 2 L 



with the proportionality constant 



8tt 



Rn 



Ln 



h 4tt i? m£ 
h c e 2 L 



(B5) 



(B6) 



depending only weakly on h/hc- By solving (jB4|) nu- 
merically and using the result in (|B6|) we can com- 
pute a, e. g. for h/hc — +1>0?— 1 its values are a = 
2.3926,2.0000,1.7838, respectively. By substituting the 
definition of A from the appendix A into (|B5[) we can re- 
late the barrier height to the experimentally known quan- 
tities Ms, d and Lo, 



E B = —fJLoM£tfL • 

47T 



(B7) 



In contrast to the collapse of a domain, nucleating a do- 
main always requires thermal activation to overcome the 
barrier formed by the exchange interaction which acts 
like a surface tension, even in absence of pinning at struc- 
tural defects. 

Splitting, or equivalently merging, of domains is hin- 
dered by an energy barrier in a similar way, since it im- 
plies a transitory increase of the domain wall length that 
is only partly compensated by a gain in dipolar energy^. 
It seems therefore plausible to assume that the height of 
the energy barrier impeding domain merger or splitting 
follows the same scaling behavior (|B7|h with a different 
numerical proportionality constant that depends on the 
exact domain configuration and the details of the pro- 
cess considered. Notice that the same scaling behavior 
E B ~ ALo has been found for the energy of a dislocation 
in a stripe pattern 3 - 3 -. 
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